Vortex nucleation as a case study of symmetry breaking in quantum systems 



G^ 
O 
O 
(N 

o 

Q 

(N 



D. Dagnino^, N. Barberan^, M. Lewenstein^'^, and J. Dalibard^ 
(1) Dept. ECM, Facultat de Fisica, Universitat de Barcelona, E-08028 Barcelona, Spain 
(2) ICFO - Institut de Ciencies Fotdniques, Pare Mediterani de la Tecnologia 08860 Barcelona, Spain 
(3) ICREA- Institucio Catalana de Recerca i Estudis Avangats, E-08010, Barcelona, Spain 
(4) Laboratoire Kastler Brossel, CNRS, UPMC, Ecole Normale Superieure, 24 rue Lhomond, 75005 Paris, France 

(Dated: December 21, 2009) 

Mean-field methods are a very powerful tool for investigating weakly interacting many-body systems in many 
branches of physics. In particular, they describe with excellent accuracy trapped Bose-Einstein condensates. A 
generic, but difficult question concerns the relation between the symmetry properties of the true many-body 
state and its mean-field approximation. Here, we address this question by considering, theoretically, vortex 
nucleation in a rotating Bose-Einstein condensate. A slow sweep of the rotation frequency changes the state of 
the system from being at rest to the one containing one vortex. Within the mean-field framework, the jump in 
symmetry occurs through a turbulent phase around a certain critical frequency. The exact many-body ground 
state at the critical frequency exhibits strong correlations and entanglement. We believe that this constitutes a 
paradigm example of symmetry breaking in - or change of the order parameter of - quantum many-body systems 
in the course of adiabatic evolution. 
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PACS numbers; 03.75.Hh, 03.75.Kk, 67.40.Vs 

In classical physics, examples of the usefulness of mean- 
field theory go back to the "molecular field theory" of mag- 
netism In the classical world, symmetry changes (or 
breaking) are driven by thermal fluctuations, and in the stan- 
dard Landau-Ginsburg scenario are associated with increase 
of classical correlations. In quantum physics, the paradigm 
example of applicability of the mean field concerns a weakly 
interacting quantum Bose gas and Bose-Einstein condensa- 
tion |i2J. The mean-field description of the gas assumes 
that its ground state ^ is approximated by a product state 
'I'(ri, . . . , rjv) — ?/'(fi) . . . ip{fM), of essencially uncoiTe- 
lated particles forming a superfluid Bose-Einstein condensate 
with order parameter ip. 

Of particular interest for quantum gases are quantum phase 
transitions and symmetry changes/breaking driven by quan- 
tum fluctuations. A celebrated example is the superfluid to 
Mott-insulator transition of bosons in an optical lattice ||3|]. 
Another example yet to be explored experimentally is the case 
of a fast rotating gas, when the number of vortices is similar 
to the number of particles, or equivalently angular momentum 
L ~ N"^ [4]. The ground state of the system is then a strongly 
correlated quantum liquid such as the Laughlin state, analo- 
gous to those emerging in quantum Hall physics [iSJ. Here, 
we consider another situation, dealing with the case of a rel- 
atively slowly rotating gas at the threshold of the nucleation 
of the first vortex. We show that owing to the symmetries of 
the system, the many-body state at nucleation is strongly cor- 
related and characterize its properties. 

The symmetry change/breaking that results from vortex nu- 
cleation has drawn a lot of attention since the discovery of 
superfluids |6J. For quantum gases, atoms are usually con- 
fined in an isotropic harmonic trap and experience an extra 
quadratic potential rotating at angular frequency Q, (for a re- 
view see ref.7). From a theoretical point of view, the vortex 
nucleation can be tackled by several techniques, ranging from 
a mean-field approach based on the Gross-Pitaevskii equa- 



tion 1181 1 10(1 to the investigation of the many-body energy 
eigenstatesllll [l2i [ll [l4l [tsj fie*. TTl . Within the mean field 
framework, standard textbooks [i2J associate vortex nucleation 
with thermodynamic instabiUty. Above a critical rotation fre- 
quency Q,c, the odd solution %jj of the Gross-Pitaevskii equa- 
tion with a single vortex [18. .191 has a lower energy than the 
even solution corresponding to the Bose-Einstein condensate 
at rest [20^. Here, we go beyond the mean-field approach and 
study the exact quantum dynamics of a mesoscopic sample of 
atoms, in the presence of the stirring potential. Our main re- 
sult is that for a rotation frequency close to Q,c, the mean-field 
description is invalid. The system enters a strongly correlated 
and entangled state, well described by an effective two-mode 
model. We compare our results with those obtained from a 
mean-field description and show that the latter exhibits dy- 
namical instability and hysteresis. As we explicitly include 
here an anisotropic stirring potential, the present mechanism 
concerns a discrete parity symmetry breaking. Therefore, it 
differs from the case of the vortex nucleation in axially sym- 
metric traps: in the latter case, breaking of the continuous ro- 
tational symmetry involves a gapless Nambu-Goldstone mode 
ll2l1l . whereas here we deal with a gapped system. 

Model. We consider a mesoscopic sample of N bosonic 
atoms of mass M placed in an axially symmetric harmonic 
potential Vb, with frequency lo±_ in the xy plane and uj^ along 
the z axis. Here, hwz is large compared with the interaction 
energy so that the dynamics along z is frozen and the gas is 
effectively two-dimensional (2D) at sufficiently low tempera- 
ture. The gas is set in rotation using an anisotropic quadratic 
potential V in the xy plane, rotating at angular frequency 57 
around the z axis. In the rotating frame, this stiiTing potential 
reads V{x, y) = 2AMu}\{x'^ — y"^), where the coefficient A 
(<C 1) measures the strength of the anisotropy. 

For A <^ 1 and 17 ~ oj^, the single-particle energy lev- 
els in the rotating frame are grouped in Landau levels, sepa- 
rated by h{Lo±^ + f7) (refs.7,22). We assume that h{uj±^ + f7) is 
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large compared with the interaction energy, so that the atomic 
dynamics is restricted to the lowest Landau level (LLL). For 
yl = 0, a basis of the LLL single-particle states is the set 
ipm{x,y) oc (a; + iy)™-e~'^^ +^ where m > is an 

integer and — yJh/Mu)i_. Each (p,,,, is an eigenstate 
of the z-component of the single-particle angular momentum 
(eigenvalue mh) and of the single-particle Hamiltonian with- 
out anisotropy (eigenvalue h[Lo±^ + m{uj± — ft)]). Within the 
LLL, we model the atomic interactions by a 2D contact po- 
tential U{f) — {h?g/M) (5(r) where g = ^/Sna/Xz is dimen- 
sionless, a is the 3D scattering length and Xz = ^/h/Mujz- 
We choose A^, hLLj± and uj± as units of length, energy and 
frequency. 

Energy spectrum. We first recall some important proper- 
ties of the A^-particle system in absence of anisotropy (A — 
0). In this case, the total angular momentum operator L com- 
mutes with the Hamiltonian H so that one can look for the 
eigenstates of H within subspaces £l of fixed L. The lowest- 
energy state in each £r. for 2 < £ < TV is lITij ITslUll : 



$L(ri, . . . ,f/v) oc ^ (Uii - Uc) . . . (Ui^ - Mc) *0 
l<ii...<ii, 

where Uj — Xj + iyj, Uc — J2j Uj/N and 

$o(?^i, . . . ,nv) oc e^^j'^j/^ . 

The energy of the state $i, is + (1 - n)L + gN{2N - 
L - 2)/(87r). At f]i = 1 - gN/iSn), all $l states for X = 
and 2 < L < iV are degenerate. The angular momentum of 
the ground state Lqs{^) shows sharp steps at critical values 
Hi, i — 1,2, .. (ref.23). Below ili, the ground state is the zero 
angular momentum state $o- At fti, Lqs jumps from to N. 
Above rJi the ground-state angular momentum has a plateau 
L = iV up to where a second jump takes place. From 
this value, a sequence of jumps and plateaux emerges up to 
the last possible L value, L = N{N — 1), corresponding to 
the Laughlin state. In the following, we focus on the vicinity 
of the first jump 57 ~ 57i, where the first vortex is nucleated. 

We now turn to the case where the rotating anisotropy 
is present. The many-body energy spectrum is calculated 
numerically by diagonalization of the Hamiltonian (see the 
Methods section). We show it in Fig.l for both zero 
anisotropy and for A = 0.03, using = 6 for illustration. 
The interaction coupling g = 1 so that 51 1 = 0.761. For 
A the ground state does not show any degeneracy around 
ill, contrary to the case A — 0. In Fig. 2, we compare LQs{fl) 
for A = and A = 0.03. For A Lqs evolves smoothly 
from to around 57 1. 

Failure of the mean-field approach for 51 ~ fJi. We now 
explain why a mean-field description must fail at 51 ~ 51 1. We 
notice that the total Hamiltonian is parity invariant. Conse- 
quently, one can look for an eigenbasis of the A^-body Hilbert 
space composed of either even or odd states. From the ground 
state of the Hamiltonian, we can extract the single-particle 
density matrix (SPDM) n'^^ (r, f) (see the Methods section). 
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FIG. 1: Energy spectrum as a function of Q,. a, Anisotropy pa- 
rameter A = Q. \), A = 0.03. In both cases, = 6 and g = 1. 
For A = Q, the ground state is multiply degenerate at the rotation 
frequency Q,\ = 1 — gN /{8n), which corresponds to the nucleation 
frequency of the first vortex. A non-zero anisotropy parameter lifts 
the degeneracy of the groun state. Here, we plot only the first nine 
energy eigenvalues from the subspace formed with even values of 
the total angular momentum, which are the only relevant ones for the 
problem addressed in this article. The arrows mark the value of Qi . 



which is also parity invariant. Hence, the single-particle or- 
bitals ipk, which are eigenstates of n^^^ with eigenvalues rife 
(J2k ''^k = N), can also be chosen with even or odd parity. 
Suppose that we vary 51 from an initial value 51^ (51; < 51i) 
to a final value 51/ (51i < 51/ < 5I2), choosing 51^./ in a re- 
gion where the mean-field description is valid, that is, when 
the largest eigenvalue ni is close to A^. For 51^ < 51i the most 
(second most) populated state V'l ("02) has no (has a) vortex 
in its central region and is even (odd). Choosing 51^ — 0.7, 
we plot the phase profiles of ^1^2 in the first row of Fig. 3 for 
A^ = 6 atoms, 5 = 1 and A = 0.03. On the other hand, at 
51 / the ground state has a single well-centered vortex and t/ji 
and ^2 are odd and even, respectively (see last row of Fig. 3 
for 51 = 0.8). Hence, the parity of ipi must change at some 
intermediate fie, which is close (for small A) to the vortex 
nucleation frequency 51i in absence of anisotropy. By conti- 
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FIG. 2: Variation of the angular momentum with rotation fre- 
quency Q. The black and red lines show, for an anisotropy A — 
and A — 0.03 respectively, the angular momentum of the ground 
state for a system of A'^ = 6 particles and an interaction strength 
g = 1. The green line is the average angular momentum predicted 
by the mean-field treatment when is linearly ramped from to 0.85 
with a slope f2 = 10"*. The initial state at f2 = is given by a slight 
perturbation of the coefficients ao = 1, ai — a2 — 0. It presents 
a dynamical instability of the zero-vortex mean-field solution for 
ft — 0.788 (marked by an arrow). Inset: The evolution of Iciip, 
which explicitly shows the instability. The blue curve of the main 
figure is the backward evolution corresponding to an initial state at 
Q — 0.85 close to the stationary mean-field solution ao — a2 = 
and ai = 1. This solution ceases to exist for Q, < 0.764, causing the 
large oscillations in the evolution of the angular momentum. 



nuity, the two most populated eigenstates tpi and -02 of n^^^^ 
must have equal populations, heralding a failure of the mean- 
field at ric- 

We show in Fig. 4 the variation of ni/N and n2/N as a 
function of n, for N = 12, g = 0.5 and A = 0.03. These 
two populations are equal for ilc — 0.775. We see that 
ni + n2 ~ N over the whole range of frequencies of this 
figure, indicating that most of the population of the SPDM is 
concentrated in the first two modes V'l and ip2- We checked up 
to iV = 20 that this concentration increases with N. Another 
relevant fact is that only the first three LLL single-particle 
states (to — 0, 1, 2) have a significant weight in the expan- 
sion of ipi and More specifically, below ilc V'l is approx- 
imately a coherent superposition of ifQ and (p2, correspond- 
ing to two off-centered vortices (even parity), whereas iIj2 is 
very close to a well-centered single-vortex state (pi (odd par- 
ity). Above ilc, tpi and ip2 abruptly exchange their form (see 
Fig.3). 

The failure of the mean-field description around ilc may 
occur in two ways. A first possibility is that for fl = ilc, the 
many -body ground level itself has a two-fold degeneracy with 
two eigenstates of opposite parity. This scenario corresponds 
to a first-order transition. It occurs when N is odd, because 
the graund state evolves from ~ Tpf^ with ipi even to ipf^ 
with tpi odd. The second possibility is that the many-body 
graund state |^'o) remains non-degenerate, as this is the case 
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FIG. 3: Density of the ground state and phase maps of tjji and 

■02. Four different values of Q for N = 6, A — 0.03 and g — 1 are 
considered. First row: Q = 0.7, ni = 5.85, n2 = 0.12. Second 
row: n = 0.760, ni = 5.01, 712 = 0.60. Third row: Q = = 
0.776 ni = 712 = 2.88. Fourth row: Q = 0.8 withni = 4.24, n2 = 
1.07. The first column is the contour plot of the total density, and the 
second and third columns show the local phase maps of tpi and tp2 
respectively. Vortices are localized at the singularities of the phase 
maps, surrounded by diffuse change of the phase. This figure shows 
that the nucleation of the first centered vortex in a rotating condensate 
by a slow frequency sweep does not occur through a smooth entrance 
of the vortex. The system passes through a correlated, non-mean- 
field state where two single-particle states have equal weight. At this 
point, -01 changes from being a coherent superposition of <^o and (p2 
(two off-centered vortices) to the single (pi state, which corresponds 
to a well-centered single vortex. Simultaneously, i/)2 experiences the 
inverse change. 



in Fig. lb. In this case, \'i'o) is even over the whole range 
[Slj, rif]. This occurs for even N and will be of interest for 
the rest of the article. 

Quantum correlations for Q. ^ Vic- We have carried out 
a detailed study of the ground state |^'o) around the critical 
frequency Vic, where the two largest eigenvalues of the SPDM 
are equal (ni = n2). At criticality, the system is very well 
described by a two-mode approximation implied by Fig.4a. 
The two largest eigenvalues of the SPDM are much larger than 
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FIG. 4: Structure of the ground state, a, Variation of the relative 
populations ni/N and n2 /N of the two most occupied states tpi and 
tp2 of the SPDM. When Q is sufficiently different from fie, ni ~ A'^, 
the system is well described by a single mode and the mean-field de- 
scription is valid. Conversely, for Q, ~ Qc, the two populations are 
comparable, corresponding to the case where a two-mode approxi- 
mation is valid even in the entangled region, b. Analysis of the state 
of the system at the critical point where ni = n2, in terms of the 
square of the scalar products Pn =\ {n : ipi ; N ~ n : tp2\'i'o) 
We obtain | {E \ *o) |= 0.92 (see equation (1)). Both panels are 
plotted for iV = 12, g = 0.5 and A = 0.03. 



all of the others, so that ni — n2 — N/2. For example, for 
N = 12, g = 0.5 and A = 0.03, we obtain ni = 712 = 
0.49 iV at fie = 0.776. The ground state is strongly correlated 
and is well described (for even N) by 



\E) = [\N, 0) + \N~ 2, 2) + ... + |0, N)]I^Nj2 + l, (1) 

where |n, m) is the state with n (respectively m) atoms in ipi 
(respectively ■4>2). Amazingly, the form of the ground state 
at VLc is practically independent of A, as long as ^4 <C 1. 
For a quantitative comparison of the exact ground state with 
the state (1), we show in Fig. 4b the squared scalar products 
{n,N — ri|5'o) in the case N = 12. They are all zero for 
odd values of n (as expected from the parity of |5'o)) and 
approximately constant for even values of n. We compared 
also our ground state at Vic with other celebrated correlated 
states, such as "Schrodinger cat" states (|iV, 0) + |0, N))/\/2, 
or "twin" states \N /2,N 12), and found much smaller over- 
laps. Although there are various ways of defining entangle- 
ment for identical particles (for a review see ll24ll ). according 
to Zanardi's concept of mode entanglement (ref.24), the state 
(1) is maximally entangled. This is clearly seen by tracing the 
state (1) over one of the two modes and observing that the von 
Neumann entropy of the reduced density matrix reaches the 
maximal value S ^ log(A^). 



At this point we mention related work on rotating ring lat- 
tices and Josephson junctions |l2^. There, strongly correlated 
states are predicted at critical rotation, but the mechanism of 
their generation, as well as their nature are fundamentally dif- 
ferent. The starting situation of these discretized models is 
that there are two degenerated single-particle states. Interac- 
tions lift the degeneracy in the many-body system and favor 
the "cat" states. In our case, the ground state for A — Q is, 
macroscopically degenerated in the presence of interactions. 
The degeneracy is lifted here by the anisotropy, leading to an- 
other kind of strongly correlated ground state. 

Vortex nucleation with adiabatic passage. We now study 
the real-time dynamics of the system using the time- 
dependent Schrodinger equation. A quasi-adiabatic evolu- 
tion that brings the system from the zero-vortex to the one- 
vortex state, is realized by sweeping il: Q,{t) = fii + 7i 
from the initial frequency Q,i chosen well below fic (typically 
Hi = 0.65) to the final frequency ilf, well above ftc (typi- 
cally VI f ~ 0.85). This evolution produces as an intermediate 
step the strongly correlated state (1). The key parameter for 
the success of this quasi-adiabatic evolution is the energy gap 
A between the ground state and the first excited state of the 
system. 

We have carried out a study of this gap for various N, keep- 
ing the product Ng constant so that Vti also remains con- 
stant. We found that for small A values (below 0.1), the gap 
is roughly constant over the range 10 < < 20, and equal 
to ~ 0.5 A. Knowing the gap, we estimate the largest pos- 
sible 7 compatible with adiabatic evolution following ref. 27 
and find 7max — £, A'^/N, where f ^ 1 (see the Methods 
section). This criterion agrees well with our results. Defining 
as successful an adiabatic evolution that leads to an overlap 
larger than 0.98 between the final state and the ground state 
at Qf, we find ^ ~ 0.1 for 10 < iV < 20. Such a quasi- 
adiabatic evolution anables us to attain the correlated state (1) 
with comparable overlap. For practical implementations, the 
atoms can be confined in a relatively tight trap at the nodes 
of an optical lattice with ujj_/2tt in the 10 kHz range. For an 
anisotropy A = 0.1 and = 10 atoms, the sweep time has to 
be of the order of one second to ensure adiabaticity. 

A natural question is the generalisation of the present 
scheme to large N. Assuming that the gap protecting the 
ground state remains constant, the mechanism will in prin- 
ciple survive. However, we have neglected here any parity- 
breaking perturbation in the Hamiltonian. Such a term would 
couple the subspaces corresponding to even and odd L val- 
ues. As shown in ref. 17, the lowest energies of these two 
subspaces are exponentially close when N increases, wich af- 
fects the robustness of the ground state. This coupling thus 
constitutes an important decoherence mechanism for large N, 
whereas our scheme remains valid for N not exceeding a few 
tens. 

Mean field approach. As our results point out that 
strongly correlated states may be reached in the course of 
the time evolution, it is interesting to see what the predic- 
tions of usual mean-field theory are. To this aim, we ex- 
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pand the condensate wavefunction f{r,t) into the relevant 
single particle LLL orbitals Lpra (^0 with angular momentum 
m = 0,1,2, /(f,i) = Em==«-?m(i)v'm(r). Using the dy- 
namical variational principle [28], we derive Lagrange equa- 
tions for the complex amplitudes a,„(t) (see the Methods 
section), and look for the stationary solutions of the form 
(im{t) — exp{—iij,t)am{0) and their stability. Finally, we 
evolve the mean-field equations and compare the results with 
the full quantum treatment. 

We choose gN = 6 and A = 0.03. Among the several pos- 
sible stationary solutions, two of them are relevant. The first 
one fa corresponds to the "no vortex" situation with a small 
admixture of "two-vortex" orbital, that is, |ao| ~ l,|a2| <C 1 
and fli — 0. This solution is the ground state for il < (l = 
0.773. The second relevant solution contains a non-zero 
contribution of the one-vortex state (ai ^ 0) and it is the 
ground state for 57 > fi. Thus, ft marks the critical value of ft 
for the thermodynamical stability of a centred vortex, and the 
"first-order transition" within the mean-field approach. 

In the frequency range 0.764 < Q < 0.788, both solu- 
tions exist and are stable, leading to a bistable and hysteresis 
behaviour. For 51 > 0.788, fa becomes dynamically unsta- 
ble (|ai I grows exponentially in time, starting from noise, see 
inset in Fig. 2). For 51 < 0.764, /{, does not exist. The nu- 
merical study confirms this hysteresis behaviour, as shown in 
Fig. 2. The green line shows the angular momentum when 51 
is ramped linearly in time from 51^ = to 51/ = 0.85, with 
the rate 51 = 10^^. A turbulent behaviour occurs once 51(t) 
reaches the edge of the stability domain of fa- The blue line 
shows the reverse evolution in which 51 varies from 51 / to 51; 
at the same rate. Evidently, the adiabatic character of the dy- 
namics cannot be maintained, in contrast to the result of the 
exact many-body treatment. 

Summary We conjecture that the scenario presented 
above is generic for the following situations: (1) it concerns 
quantum mechanical systems in which the ground state un- 
dergoes symmetry change/breaking as some parameter of the 
system A crosses a critical value Ac; (2) far from Ac, the sys- 
tems are well described by the mean-field theory with order 
parameters reflecting the change of symmetry; (3) in the dy- 
namical mean-field description, the system exhibits dynami- 
cal instability and breakdown of adiabaticity. 

In such situations we expect the appearance of strongly 
correlated states. The SPDM shows typically a few rele- 
vant single-particle modes that are involved in the symmetry 
change. They can be guessed by analyzing the results of the 
dynamical mean-field approach. For instance, if this approach 
exhibits standard signatures of bistability, we can expect two 
relevant modes as in the case study presented here. Similar 
insight can be gained from analysis of small Gaussian fluctu- 
ations around the mean-field solutions, that is, Bogoliubov- 
de Gennes equations j29ll. Reduction of the full theory to 
the quantum modes provides then a very good approxima- 
tion. Alternatively, it can be viewed as re-quantization of the 
mean-field theory reduced to the relevant single-particle or- 



bitals (vf\. The strongly correlated states appearing in such a 
situation exhibit strong entanglement and this property can be 
detected in experiments with moderate N. 

Methods 

Diagonalization of the Hamiltonian. In the frame rotat- 
ing at angular frequency 51, the Hamiltonian of the system is 
H = Hq + U, where Hq is the sum of one-body Hamilto- 
nians Hq — J2f=i ^o.j and U is the two-body interaction 
potential, characterized by the 3D scattering length a. Each 
one-body Hamiltonian is the sum of kinetic, potential and ro- 
tation energy: 

uj 2M ' 2 ^ ^ ^ II ■'^ 

- 51^2^ + Vj 

where Vj — 2 AM uj\ (a;| — j/| ) is the anisotropic potential that 
sets the gas in rotation. We assume that the interaction energy 
is much smaller than hwz so that the z motion is frozen and 
the atoms occupy only the ground state exp(— 2:^/(2A2)) of 
this degree of freedom. The gas is supposed to be rotating 
sufficiently fast to have wj, — 51 <C wj. + 51, which guaran- 
tees that the various Landau levels are well separated from 
each other. The interaction energy is also assumed to be small 
compared to K(lij^ + 51) so that the low temperature dynamics 
is restricted to the LLL. 

In the absence of anisotropic potential A — the eigen- 
states of the one-body Hamiltonian in the LLL are the func- 
tions (p,n{x,y) (X {x + iy)'"e-(^'+2'')/2^i, m = 0, 1,2, . . .. 
We introduce the creation a}^^ and annihilation operators 
of an atom in state (p„i, and we write H in the second quanti- 
zation 

H = h uj^N + h {LU^-n)L + V + U , 

where N ~ J2o-m'^m and L = ^TOaJ,ja„i are the particle 
number operator and the total z-component angular momen- 
tum operator, respectively. The expression of the rotating po- 
tential in the second quantization is 

m 

+ V(to + 1)(to + 2) al, • 
Finally the contact interaction potential reads 

where the matrix elements are given by 

U1234 = (toi m2 I [/ I mg 7714) 

_ g I^mi+m2,m3+r«4 (mi + 7712)! 
~ A^TT VTOl!m2 17773 17774! 2™l+"'2 + l 
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In the absence of anisotropy (A = 0), H and L com- 
mute and share a common basis. The first step in the diag- 
onaHzation of the Hamihonian is to determine a basis |Ap) 
(p = 1, . . . , Tii) for each subspace of given total angular mo- 
mentum L. The dimension of each subspace corresponds 
to all of the possible configurations of N particles with an- 
gular momentum rrij that fulfil the condition L = EjLi "^j- 
The matrix of the Hamiltonian in the LLL basis then consists 
of blocks of size x n^, which we diagonalize using stan- 
dard codes. 

When A 0, the anisotropic potential connects the various 
subspaces of given L. We then choose a maximum angular 
momentum imax and write the matrix giving the restriction 
of the Hamiltonian to the subspace of states with L < L-^^^. 
This Q X Q matrix, with Q = J2l=o° again diagonal- 

ized using standard codes. In practice the value of L,„ax is 
chosen to ensure a good convergence for the energies and the 
eigenstates of the Hamiltonian. The results given here have 
been obtained with L„iax = + 2. 

Note that the anisotropic rotating contribution V can in 
principle be included within the framework of the LLL ap- 
proximation in two ways. The first approach has just been de- 
scribed above and consists of keeping the same Landau levels 
as for A = and then diagonalizing H within the LLL. The 
second approach consists of calculating exactly the single- 
particle eigenstates in presence of the anisotropy V, and defin- 
ing a new LLL accordingly llsoll . The Hamiltonian is then 
diagonalized within this 'anisotropic' LLL. We have checked 
that both methods lead to very similar results for $7 ~ 
The results presented here have been obtained with the first 
approach. 

Single particle density matrix The SPDM can be re- 
garded as an integral operator with the kernel: 

n^^\f,r') = {^a \ *^(f) #(r')|*o), 

with ^(r) and being the annihilation and creation field 
operators of an atom in r. The single-particle orbitals are the 
eigenstates of the SPDM: 



If there exist a single relevant eigenvalue such that ni ^ 
X]fc>2 "-fc' ^^^^ y^i^ii^ has the role of the order parame- 
ter of the system. In particular, the map of the local phase 
of this function gives precise information on the location of 
vortices 1 15]. 

Adiabatic approximation The diagonalization of the 
many-body Hamiltonian provides the eigenstates |5'j (r2)) and 
the eigenenergies Ej{il,). In particular, the ground state 
|5'o(^l)) is separated from the first excited state by an 

energy gap Hluiq (fl), which is minimal at the avoided crossing 
close to ill. We consider here a process where is scanned 
linearly from il^ < to fi/ > fii and we want to find a 
criterion on fl ensuring that the system follows adiabatically 
the ground state, with a negligible transition rate to the other 
states. 



The probability for a non-adiabatic transition is 
given by It27i1 : 



Po^j 5: rnax 



where ajo = (*j|(d|^'o)/dt). We have 



rfl^o) 
dt 



dn 



From the eigenvalue equation H\"^q) = i^ol^'o). we obtain 
after a derivative with respect to fl: 



d\^o) dEo 



dn 



dn 



*o(f^)) + Eo 



dfl 



We now take the scalar product with (^fj | {j ^ 0) and we get: 



dn 



We choose equal to the first excited state of the system 
l^&i). The matrix element (^I^i |Lz|\I>o) is at most of order Nh 
in the vicinity of the avoided crossing. Therefore: 



aio 



c#o> 
dt 



< n 



Nh 



fkui 



hence the condition for po^i ^ 1: 



■ N 

O ^ < 1 



Mean-field approach. The mean-field approach consists 
of assuming that all atoms are in the same state f{f,t) = 
Em=o a.ni{t)ipm{r) with J2 lomP = 1- The average angular 
momentum per particle is L — \ai\ 



2 1 a2 1 and the average 
energy per particle E{Tp) — j^{f'^^\H\f'^'^) reads (up to an 
additive constant): 

E{ij) = {l-n){\ai\^ + 2\a2\^) + V2A{aoa*2+a*a2) 



+ 



Ng 
47r 



"2 



+2|aonai|' + |aona2|' + Tr|aina2| 



'V2 



The Lagrange equations associated with this energy are idj 
dE/da* (ref. 28), which gives for example: 



^PiA 



Ng 
27 



ao(|ao|^ + |ai|^ + ^|a2|^) + 



1 



2\/2 



2 * 
^10-2 



and two similar equations for di and a2- Note that in this 
mean-field approach, N and g have a role only through the 
product Ng. In particular, the fact that N is even or odd is of 
no relevance here. 
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The stationary solutions are obtained by inserting a„i{t) = 
am(0) e"*''* in the three Lagrange equations. A detailed anal- 
ysis of the resulting 3x3 nonlinear system shows that two 
classes of solution exist. The first class (/„) corresponds to 
ai = 0. Depending on the value of the parameters Ng,A and 
Cl, there may exist two, three or four solutions of this kind. 
After some tedious but straightforward calculation, one can 
obtain for this first class of solution an analytical relation be- 
tween fl and the angular momentum per particle i = 2|a2|^: 

Stt V 8 ; ^L{2 - L) 

The second class of solution corresponds to a non-zero value 
for a\ and we have not been able to provide an exact analytical 
expression for the solution in this case. Using a numerical 
analysis, we have determined the local minima of the energy 
and we found that one solution of this kind exists if and only 
if f2 > 0.766. We have compared the energy of this solution 
with the lowest energy of the solutions in the first class: for 
12 < = 0.773 (respectively Q > fi), the the ground state 
is obtained with a solution belonging to the first (respectively 
second) class. 

The stability of the solutions of the first class (ai = 0) can 
be studied analytically by looking at the equation of evolution 
of &i = aie*''*. This equation can be linearized around 6i = 
and written in the form ihx = Abi + Bb^, where the constants 
A and B are real numbers that can be calculated explicitly in 
terms of the parameters li, A and Ng. The stationary solution 
corresponds to bi = 0, and it is stable if bi {t) stays around 
when starting from a small non-zero initial value. This hap- 
pens when 1^1 > \B\, whereas bi undergoes an exponential 
divergence from any initial noise ii \A\ < \B\, signalling a 
dynamical instabiUty of the solution. 
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